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ABSTRACT 

Context. Intermediate mass protostars provide a bridge between low- and high-mass protostars. Furthermore, they are an important 

component of the UV interstellar radiation field. Despite their relevance, little is known about their formation process. 

Aims. We present a systematic study of the physical structure of five intermediate mass, candidate Class protostars. Our two goals 

are to shed light on the first phase of intermediate mass star formation and to compare these protostars with low- and high-mass 

sources. 

Methods. We derived the dust and gas temperature and density profiles of the sample. We analysed all existing continuum data on 

each source and modelled the resulting SED with the ID radiative transfer code DUSTY. The gas temperature was then predicted by 

means of a modified version of the code CHT96. 

Results. We found that the density profiles of five out of six studied intermediate mass envelopes are consistent with the predictions 

of the "inside-out" collapse theory. We compared several physical parameters, like the power law index of the density profile, the size, 

the mass, the average density, the density at 1000 AU and the density at 10 K of the envelopes of low-, intermediate, and high-mass 

protostars. When considering these various physical parameters, the transition between the three groups appears smooth, suggesting 

that the formation processes and triggers do not substantially differ. 

Key words, stars: formation 



1. Introduction 

Stars are usually catalogued in three groups according to their 
mass: low-mass (<2 M ), intermediate mass (IM) and high- 
mass (>8 M ). The high-mass stars are usually found to form 
in tight clusters (e.g. Hillenbr and & Hartma nn 1998), including 
many low-mass stars as well. Low-mass stars on their own can 
be formed i n isolation or in loose groups of a few objects per cu- 
bic parsec (Go mez et al.lll993l) . While there is reasonable con- 
sensus on the overall process leading to the fo rmation of iso- 
lated, low-mass stars (see e.g. Shu & Adams 1987), the situation 
for high-mass star formation is much more confusing. The major 
problem is that a central object with M > 8 M© produces a ra- 
diative pressure stronger than the gravitational pressure, which 
should stop the accretion. Therefore stars with a mass greater 
than ~8 M© cannot be formed in a simple-minded fully gravi- 
tationally driven process and other mechanisms need to be con- 
sidered. Studying the differences in observed environments, sev- 
eral authors have proposed that high-mass stars are formed by 
the coalescence of lower-mass stars, while other authors favor 
"monolithic" formation (see for example the recent review by 



iBeuther et alJ 120071). where turbulence or winds/outflows c an 
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play a major role dMcKee & Tanll2003l: iKrumholz et al.ll2005l) . 

The differences between low- and high-mass protostars are 
not only restricted to the dynamics of the formation mecha- 
nism. Notabledifferences also exist when it comes to the chem- 
istry. Bottinelli (2006) compared the abundances of O-bearing 
complex organic molecules (COMs) present in Hot Cores sur- 
rounding high-mass protostars and in Hot Corinos around low- 
mass protostars. They concluded that the abundances of the 
O-bearing COMs, compared to methanol or formaldehyde, are 
about one order of magnitude higher in Hot Corinos than 
Hot Cores. In addition, significant differences in the abun- 
dances of molecules constituting the grain mantles (H2O, CO, 
C02,CH30H, O CN~) were observed between low- and high- 
mass protostars dBoogert et alJ |2004). Finally, in contrast to 
high-mass star forming regions, low-mass protostars present an 
extremely enhanc ed molecular deutera tion, the so-called Super- 
Deuteration (e.g. ICeccarelli et all 2007), not observed in high- 
mass star forming regions. 

The IM stars, namely stars whose mass lies in the 2 
to 8 M range, are crucial in studies of star formation be- 
cause they provide the link between low- and high-mass stars 
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(iDi Francesco et al.lll997l:lMannings & Sargentll997Ll2000h and 
therefore can help to determine whether different formation pro- 
cesses are required. IM stars are also found in clusters (e.g. 
iTesti & Sargenil 19981: iNeri et al.ll2007l: JFuente et alJl2007b . with 
a smooth transition towards the low- mass star, loose cluster 
regime for star masses around 3.5 M© (ITesti etalJI 19991) . These 
authors also conclude that IM stars mark the transition from 
the low-density aggregates of < 10 stars per cubic parsec of T 
Tauri stars to the dense clusters of > 10 3 stars per cubic par- 
sec associated with massive stars. Finally, IM stars are among 
the most dominant sources for the interstellar FUV field (e.g. 
lHabingll968l:lGondhalekar & Wilson! 19751) . which regulates the 
phases of the ISM in the Galaxy and, in turn, the overall Galaxy 
star formation process and history. In this context, the study of 
the formation of the IM stars is not only interesting in itself, but 
it can also greatly help the debate on the formation of high-mass 
stars. Indeed, due to their intermediate position, the study of IM 
protostars provides crucial information about the transition be- 
tween the two formation regimes, as well as on the limits of the 
low-mass and high-mass formation scenarios. 

Despite their importance, very little is known about the 
formation and first evolutionary stages of IM stars. Although 
some individual IM protostars hav e been studied in detail 
dFuente et al] 120071: iNeri et al] 120071: iHogerheiide et al] 1 19991: 
Sch reyer et al.ll2002l) . the situation is so bad that to date we nei- 
ther have a satisfying sample of Class IM objects, namely 
objects representing the first phases of stellar formation where 
the protostar is embedded in its envelope and its luminosity is 
dominated by the accretion luminosity, nor a systematic study 
of their physical structure, as in the case for low-mass and high- 
mass protostars. On the contrary, several previous studies have 
considered samples of low- and high-mass protostars. Relevant 
similar studies to the pres ent IM one are tho se reported by 
iJorgensen eFaD (120021) and IShirlev et all (120021) . who targeted 
7 and 11 low-mas s Class protostars respectively, and the sys- 
tematic studies by Mueller et al. (2002), Hatchell & van der Tak 
(120031) . IWilliams et al] (120051) and IVan der Tak et al] (12000). 
who targeted 31, 10, 32, and 9 high-mass protostars respectively. 
These authors derived the physical structure of the envelope for 
each protostar using continuum observations coupled with a ra- 
diative transfer code assuming a ID spherical model. Combined, 
these studies cover a luminosity range of 1 < L < 22, and lxlO 3 
< L < 2xl0 6 for low- and high-mass protostars respectively, 
leaving a significant gap in luminosity and knowledge. 

This paper aims to fill this gap, providing a systematic study 
of the physical structure of the envelopes of a sample of five IM 
protostars. In addition, this study will provide a valuable input 
for future work on the c hemical structure of IM protostars (e.g. 
lAlonso-Albi et al .1120091 in prep). The pre sent work follows the 
method presented in lCrimier et al] J2009a|) to study the IM pro- 
tostars OMC2-FIR4. Namely, the combination of sub-millimeter 
to infrared observations with radiative transfer codes are used 
to derive the physical structure of each envelope. We selected 
a sample of five representative candidate Class IM protostars: 
CB3-mm, Cep E-mm, IC1396 N BIMA 2, NGC7129 FIRS 2, 
and Serpens FIRS 1 . These sources were selected using the fol- 
lowing criteria: 



- the luminosity is in the 30 to 1000 L range (corresponding 
to a ~ 3 to 9 M Q of ZAMS stars); 

- the continuum maps at submillimetre and/or millimetre 
wavelengths are available; 



- the source is closer than 2.5 kpc to the Sun and isolated, 
namely no source of similar brightness in the submillimetre 
within an angular distance of -100". 

The coordinates, the distances and the bolometric luminosities 
of the selected sources are listed in Table [T] 

This article is organised as follows. Section [2] describes the 
observations used and the method and tools employed to derive 
the dust and gas density and temperature profiles. Section [3] first 
summarises the general results obtained toward the IM proto- 
star sample. Next, for each source, background information, data 
available, and the results obtained are presented. In Sect, ffl the 
physical structure of the envelopes of IM protostars obtained in 
this study are compared with those obtained in previous stud- 
ies of low- and high-mass protostar envelopes. In this section 
we also discuss the limitations of our modeling method. Finally, 
Sect. \5\ concludes the article. 

2. General method 

In this section we present the general method used to derive the 
dust density and temperature profiles by modelling the contin- 
uum maps of each region, plus the Spectral Energy Distribution 
(SED) from the millimetre to the Mid-Infrared (MIR). We first 
describe the observations used in the analysis and how the 
brightness profiles and the SED are derived (Sect. 12. ll) . and then 
we present the modelling (Sect. 12.21) . 

2.1. Continuum emission : observational data 

In our analysis we used the maps of the continuum emission 
at 850 and 450 pm obtained at the 15 m James Clerk Maxwell 
Telescope (JCMT). In addition, we took into account the SED 
of each source from 24 /mi to 1300 /mi data obtained with the 
Spitzer, IRAS, JCMT and IRAM telescopes. Note that specific 
details (datasets, pipeline, uncertainties, contamination ... etc) 
are given separately for each source in Sect. [3j 
a ) Continuum brightness profiles 

For each source we retrieved from the JCMT Science Archivqj 
the 450 yum and 850 /mi maps obtained by the focal-plane instru- 
ment SCUBA (Submillimetre Common-User Bolometer Array). 
The data were reduced with the JCMT pipelines SURF VI. 6- 13, 
and KAPPA 1.9-7. When multiple observations were available 
in the archive, we used the most sensitive ones. The angular res- 
olution of the maps is approximately 7.5"and 14.8"at 450 and 
850 pm respectively (see Sect. 12.21 for more details). Based on 
the many previous JCMT published observations, the calibra- 
tion uncertainty and noise levels are < 10% and 0.04 Jy beam -1 
at 850 yum and < 30% and 0.3 Jy beam -1 at 450 pm, respectively. 
To evaluate the continuum brightness profile (the brightness as 
a function of the distance from the envelope centre), we aver- 
aged the observed flux over annuli with the same distance from 
the centre of the source. The final uncertainties were evaluated 
taking into account the calibration uncertainty, the noise levels 
of the used maps, and the non-sphericity of the source. The non- 
sphericity of the source is taken into account via the computa- 
tion of the standard deviation relative to the azimuthal average 
for each averaged flux of the brightness profile. 
b)SED 

The SED points at 850 /mi and 450 /mi were obtained integrating 
the continuum flux profiles previously derived over the full enve- 
lope. We also considered the integrated flux at 24 pm and 10 pm 



1 http://www.jach. hawaii. edu/ JCMT /archive/ 
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extracted from the Spitzer Space Telescope Multiband Imaging 
Photometer (MIPS) maps. To this end, we retrieved the obser- 
vations from the Spitzer Science archival. The data reduction 
was performed using the pipeline SI 6. 1.0. The evaluation of the 
fluxes was done with the same method employed for the maps 
at 850 jim and 450 yum, namely deriving the brightness profiles 
from the maps and integrating over each annuli. We evaluated 
the uncertainties according to the instrument uncertainties and 
varying the outer radius of the envelope to account for the uncer- 
tainty in the envelope size. To do this we first set the outer radius 
equal to where the flux profile becomes flat (because dominated 
by the noise or the flux of the parent cloud) and varied this value 
up to 100 %. In general, the uncertainties remain relatively con- 
stant for variations larger than 30-50%. Note that we checked a 
posteriori that our adopted radius and variation range is consti- 
tent with the outer radius obtained after the modelling. Finally, 
we also considered the IRAS fluxes at 60 /mi and 100 //m, re- 
spectively, taken from the literature or the IRAS Point Source 
Catalog v2.C0. The SED fluxes of each source and their uncer- 
tainties are reported in Table [2] 

2.2. Continuum emission : modelling 

To derive the dust temperature and density profiles of the 
envelope, we used th e ID radiative transfer code DUSTY 
(llvezic & Elitzur|[l997l) . Briefly, giving as input the temperature 
and size of the central object and a dust density profile, DUSTY 
self-consistently computes the dust temperature profile and the 
dust emission. The comparison between the computed 450 fim, 
850 /mi and 1300 fim brightness profiles (namely the brightness 
versus the distance from the centre of the envelope) and inte- 
grated SED against the observed brightness profiles and inte- 
grated SED (see previous paragraph) allows one to constrain the 
density profile and, consequently, the temperature profile of the 
envelope. 

To be compared against the observations, the theoretical 
emission must be convolved with the beam pattern of each tele- 
scope. Following the recommendations for the JCMT, the beam 
is assumed to be a combination of three Gaussian curves: at 850 
fim we used HPBWs of 14.5", 60", and 120", with amplitudes 
of 0.976, 0.022, and 0.002 respectively; at 450 //m the HPBWs 
are 8", 30", and 12 0". with amplitude ratios of 0.934, 0.06, and 
0.006, respectively (ISandell & Weintraubll200H) . 

In all sources we assumed that the envelope density follows 
a single index power law 



2004 lYounget aDl2003l), we adopted the dust opacity calcu- 
lated by Ossenkopf & Henning (1994), specifically their OH5 
dust model, which refers to grains coated by ice. We also ob- 
tained a lower limit to T in (the temperature at R/„) of 300 K: 
any higher value for T in would give similar results. Again, given 
the high optical thickness of the envelope at short wavelengths 
and the low contribution from the dust at T > 300 K to the sub- 
millimetre emission ( as the NIR emission is mainly from the 
very inner part, i.e. a small volume), increasing T in does not 
make a difference to the resulting best-fit. 

In summary, the output of DUSTY depends on four free pa- 
rameters: a, no, r in and x out . In practice, the DUSTY input pa- 
rameters are the power law index, a, the optical thickness at 
100 yum, noo, the ratio between the inner and outer radius, Y 
(=r out lTi n ) and the temperature at the inner radius T in . The opti- 
cal thickness is, in turn, proportional to the dust column density 
which depends on no and the physical thickness of the envelope. 
Note that because the beam sizes of the available maps are rel- 
atively large (> 7.5", which corresponds to a radius of > 3250 
AU for a source at a distance of 1000 pc), the inner regions of 
the envelopes are relatively unconstrained by the available ob- 
servational data. Finally, as explained in Ivezic & Elitzuri (1 19971) . 
DUSTY gives scaleless results (which make it very powerful be- 
cause the same grid of models can be applied to many different 
sources). To compare the DUSTY output with actual observa- 
tions it is necessary to scale the output by the source bolometric 
luminosity, L^/, and the distance to the source. Note that the 
bolometric luminosity is estimated by integrating the emission 
across the full spectrum. By definition, this can only be done 
when the entire SED is known. This is exactly one of the out- 
puts of the modelling. Therefore, we re-evaluate the luminosity 
of each source iteratively from the best-fit model, by minimising 
the* 2 . 

We ran a grid of models to cover the parameter space as re- 
ported in Tabled The best-fit model is found by minimising the 



Parameter 


Range 


step 


a 


0.2-2.5 


0.1 


Y 


50-2000 


10 


T"l00 


0.1-10. 


0.1 


^ 


300 K 




T* 


5000 K 





n(r) = no x 



w- 



(1) 



Table 3. Range of the input parameters to DUSTY covered in 
the present study. The range of the a, Y and T100 parameters are 
covered by increasing by a constant value at each step of the grid 
reported in the third column. Note that T in and T + remain fixed 
as they do not influence the results (see text). 



where the power law index, a, and the density, no, at ro are 
free parameters of the model. The envelope is assumed to start 
at a radius r in and extends to r out . Both r in and v out are addi- 
tional free parameters of the model. The last input into DUSTY 
is the temperature of the central source, T*, here assumed to 
be 5000 K. We verified that the choice of this last parameter 
does not influence the results. Indeed, given the high optical 
thickness of the envelopes at the wavelengths where the central 
source emits, the model outputs are not sensitive to the T* value. 
Finally, the opacity of the dust as function of the wavelength 
is anothe r parameter of DUSTY. Following numerous previous 
studies (Ivan der Tak et al.ll 19991: [Evans et alJkOOlHShirley et al] 



2 http://ssc. spitzer. caltech. edu/archanaly/ 

3 http://irsa. ipac. caltech. edu/cgi-bin/Gator/nph-dd?catalog-iraspsc 



X 1 with an iterated two-steps procedure. First, we used the ob- 
served brightness profiles at 450 \im and 850 \im to constrain 
a and Y, assuming a value for noo. The x^ e d computed during 
this first step are reported as ^^-maps in Table |4] Second, we 
constrain the optical thickness noo by comparing the computed 
and observed SED, assuming the a and Y of the previous step. 
The new noo is used for the next iteration and the process is re- 
peared. The x 2 re d computed during this second step are reported 
as^^-SED in Tableffl In practice, the iteration converges in two 
steps. This occurs because the normalised brightness profiles are 
only very weakly dependent on noo, while they are very depen- 
dent on the sizes of the envelope and on the slope of the den- 
sity profile (see also iJdrgensen et al J 120021: ISchoier et aDl2002t 
Cri mier et al.ll2009ab . On the contrary, the optical thickness de- 
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pends mostly on the absolute column density of the envelope, 
which is constrained by the SED. 



2.3. Gas temperature profile : model description 



CeccarellietaD d 1996b . 
Dotv & Neufeldl d 1997b s 



ICrimier et all (l2009ah 
showed that the gas can be 



and 
can be ther- 
mally decoupled from the dust in the inner regions of low-, 
intermediate and high-mass protostellar envelopes. The de- 
coupling occurs mainly in the inner part of the envelopes: for 
example whereas T dust ~ 150 K at 200 AU, T gas ~ 90 K. Note 
that in the outer envelope, where Tdust < 100K, gas and dust 
temperature only di ffer by a few percents (see discussion in 
ICrimier et al.ll2009ah . The reason for this decoupling is the high 
water abundance in the gas phase caused by the sublimation 
of grain mantles. We therefore explicitly computed the gas 
temperature profile of the envelope surrounding each source 
by finding the equilibrium temperature obtained by equating 
the gas cooling and heating t erms at ea c h radi us. Following the 
method described by Ceccare fli et al.l (Il996b . we considered 
heating from the gas compression (due to the collapse), dust-gas 
collisions, photo-pumping of H2O and CO molecules by the 
IR photons emitted by the warm dust close to the centre, and 
cosmic rays ionisation which is a minor heating term in proto- 
stellar envelopes. The cooling is mainly due to the rotational 
lines from H2O and CO, plus the fine structure lines from O. 
The gas temperature therefore depends on the abundance of 
these three species. In practice, only the water abundance is an 
important parameter of the model, because the CO and O lines 
are optically thick and LTE populated in the range of CO and 
O abundances typical of protostellar envelopes, while the water 
levels are sub-thermally populated (non-LTE). In the non-LTE 
regime the water levels are excited by collisional processes 
and de-excited by radiation while in LTE regime the levels 
are mainly excited and de-excited by collisions. Therefore, 
a photon emitted in LTE regime in an optically thick region 
will be absorbed and then the absorber will be de-excited by 
collisions while in non-LTE regime the absorbed photon will 
possibly be re-emitted. In this case we say that the water lines 
are effectively optically thin. For this reason we computed 
two cases for the water abundance, as it is generally poorly 
constrained in protostellar envelopes and totally unconstrained 
in intermediate mass protostars. We adopted a step function for 
the water abundance profile to simulate the jump caused by ice 
sublimation. The jump is assumed to occur at 100 K. We con- 
sidered the H2O abundance (with respect to H2) X(H20)/ n in the 
inner envelope, where T> 100 K, equal to 10" 5 and 10" 6 , fixing 
the water abundance in the outer region, X(H20) 0M? , at 10" 7 . 
The CO and O abundances were fixed a t the standard values 
found in mol ecular clouds, i.e . 1 x 10~ 4 ([ Frerking et al.lll982h 
and 5 x 10" 4 dCaux et al.lfl999l : I Vastel et alJl2000h . respectively. 
Note that because the O and CO lines are mostly optically 
thick the exact value of their abundance is not important. To 
c ompute the cooling from the lines we used the code described 
in lCeccarelli etafl d2003Lll996h and lParise et all d2005b . Briefly, 
the line cooling is computed with an escape probability method, 
which takes into account the dust level pumping and the line 
optical depths at each point of the envelope by integrating over 
the solid angle. A re cent description of the code is reported in 
ICrimier et al. ( 2009a. b). The same code has been used in several 
past studies, whose results have been substantially confirmed 
by other groups ( e.g. the analysis on IRAS 16293-2422 by 
ISchoier et al.ll2002h . 



For the collisional coefficients of water w it h hyd rogen 
molecul es, we used the data by iFaure et all d2007b and 
Fau re & Josselinl (2008) available for the temperature range 20- 
5000K. Because the ortho-to-para conversion process of H2 is 
chemical rather than radiative, the Ortho-to-Para Ratio (OPR) 
H2, which the water population depends on, is highly uncer- 
tain. The recent analysis of H2CO o bservations towards a cold 
molecular cloud by Troscomp t et al] d2Q09b confirms theoretical 
estimates (e.g. iFlower et al.l l2006) that in molecular clouds the 
H2 OPR is lower than 1 . Lacking specific observations towards 
protostars, here we assume that H2 OPR is in Local Thermal 
Equilibrium and therefore follows the Boltzmann distribution 



OPR 



(2/ + l)£(2/+l)exp(-^) 
(2I P + 1)2(2/ + l)exp(-^)' 



(2) 



where I and I p are the total nuclear spin, corresponding to 
whether the hydrogen nuclear spins are parallel (I = 1, TT) or 
anti-parallel (I p = 0, |J,). The sum in the numerator and de- 
nominator extends over all ortho and para levels J, respectively. 
Similarly to H2, water comes in the ortho and para forms. In 
these cases, because the water is the dominant gas coolant only 
in the regions where the dust temperature exceeds 100 K, we 
assumed an OPR equal to 3, strictly valid for gas temperatures 
higher than 60 K. Because the water lines are optically thick, 
the cooling depends on the velocity field, assumed to be that 
of an envelope collapsing in free-fall towards a central object 
with a mass M*. Here, we assumed that the entire envelope is 
in infall. The masses M* used for each source are reported in 
T able ffl They were d erived analytically, following the equation 
in Stah ler et all (1986) which links the mass of the central object 
to the bolometric luminosity, the accretion rate, the mass and the 
radius of the envelope. Basically, the equation assumes that the 
luminosity is entirely due to the gravitational energy released in 
the collapse and uses the computation of the hydrostatic core 
radius of the protostar by IStahler et al.l dl986b . We checked the 
influence of our results against this assumption, running cases 
with M* fixed at 2 M©. The difference in the gas temperature 
between the two cases never exceeds 1%. 



2.4. Water line observations 

In order to constrain the water abundance, which is very impor- 
tant for computing the gas temperature, we considered obser- 
vations obtained by the Long Wavelength Spectrometer (LWS) 
aboard the Infrared Space Observatory (ISO) in the 45 //m to 
200 yum range, where several water lines emit. All sources ex- 
cepted CB3 were observed with the LWS. In two sources, Cep 
E-mm and Serpens FIRS 1 , several water lines were detected and 
thei r analysis has been reported by More-Ma rtin et al.l d200lb 
and lLarsson et a l. (2002), respectively. As discussed by these au- 
thors, given the relatively large beam of the LWS, the water line 
fluxes are due to the combination of many components along the 
line of sight: outflows, multiple sources and Photo-Dissociation 
Regions. The measured H2O line fluxes are thus upper limits to 
the fluxes from the envelopes and we checked that our predic- 
tions do not exceed the observed fluxes. For two other sources, 
IC1396 N BIMA 2 and NGC7129 FIRS 2, we retrieved the LWS 
grating spectra (spectral resolution Av/v -200) from the ISO 
Data Archive! and extracted the upper limits to the flux of the 
brightest H2O lines. Table \5\ summarises the water line ISO ob- 
servations for each source. 



http://iso. esac. esa. int/ida/ 
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3. Results : Dust and gas density and temperature 
profiles 

The general results of our analysis are: 

- The continuum brightness profiles and the SEDs at wave- 
lengths larger than 60 jim for the five IM-proto stars of our 
sample can be reproduced by spherical, single index, power 
law density models. On the contrary, in all sample sources, 
this class of models fails to reproduce the 24 /mi flux, un- 
derestimating it by 1 to 3 orders of magnitude. The possible 
causes and implications of this failure are discussed in detail 
in Sect, ffl The ^^-maps and ^^-SED obtained for each 
source are reported in Table ffl The ^^-maps ranges from 
0.1 1 to 0.55 for the sources modelled in this paper. Th&x^ ed - 
SED ranges from 0.8 to 4.4 for the sources modeled in this 
paper. Note that because the flux at 24 yum is underestimated 
by several orders of magnitude, the^^-SED value is mainly 
driven by this point. 

- The power law index for the five sample sources varies be- 
tween 1.2 and 2.2, with an average value equal to 1.6. 

- The envelope radius varies between 6000 AU for the lowest 
luminosity and closest source, Serpens FIRS 1, and 1 x 10 5 
AU for the brightest and farthest source, CB3-mm. 

- The radius where T^ ? =100 K lies between 100 (Serpens 
FIRS 1) and 700 (CB3-mm) AU. 

- The density at T^=100 K varies, from 0.4 to 2 x 10 8 cm -3 . 

- The envelope mass ranges from 5 M© (Serpens FIRS 1) to 
120 (CB3-mm) M . 

- The mass of the central object is estimated to be between 0. 1 
M (Serpens FIRS 1) and 6 (CB3-mm) M Q . 

- The gas and the dust are thermally coupled across the en- 
velope in three of our sources, with temperature differences 
lower than 5 K. For IC1396 N BIMA 2 (Sect. E3 and 
NGC7129 FIRS 2 (Sect. El) the difference can reache 40 
K. 

- The predicted H2O lines are consistent with the ISO upper 
limits of Sect. [231 



Table H summarises the best-fit parameters and some relevant 
physical quantities derived from the dust radiative transfer anal- 
ysis of each source, Table[5]lists the water line predictions. Note 
that given the relatively large beam, the ISO observations are 
contaminated by the emission from outflows, multiples sources, 
and Photo-Dissociation Regions and are therefore considered 
only as upper limits. 

The appendix describes in detail the source background, the 
data included in the analysis, and the derived physical structure 
(gas and dust density and temperature profiles) for each source. 



4. Discussion 

4.1. The link between low- and high-mass protostars 

One of the major goals of this work is to verify whether interme- 
diate mass protostars provide a link between low- and high-mass 
star formation. In this section, we analyse whether the parame- 
ters describing the protostellar envelope structure (power law in- 
dex, dust temperature at a given distance, envelope mass...) de- 
pend on the luminosity, and hence the mass, of the future central 
star. 

Figure [T] plots key parameters of the envelope structure 
(the power law index of the density profile, a, the total mass, 
M env , the outer radius, r ouU and the average density, <n> 



of the envelope) for low-, intermediate, and high-mass pro- 
tostars as a function on the bolometric luminosity of each 
source. The ensemble of the plotted sources covers six orders 
of magnitude in luminosity, from about 1 L© to ~ 10 6 L . 
The plotted d ata are from the present study (Table |4]) and 



Crimier et al.l (2009a) for the intermediat e ma ss protostars, 
Jorgense n et al.l (120021) . IShirlev et al] (120021) and ICrimier et al. 
(2009b) f or the low-mass proto stars, and IVan der Tak et al. 



(2005) 



(2000), Hatchell & van der Tal3 (120031) . IWilliams et al 
|5|) and Mueller et al.l (l2002h for the high-mass proto - 



stars. Note that lJ0rgensen et al.l (120021), [W illiams et~all (12005 

"TF" 



lHatchell & va n der Tak (2003), and Crimier etal. (2009a.b), 
as well as the present study, use the DUSTY code in the 
analysis. The correlation coefficients and probability for a 
chance correlation between the pair of parameters considered in 
Fig. [T] are reported in Table The plots and quantities shown 
in Fig. [T]and Tabled respectively, lead to the following remarks. 

Density power law index a: 

The density power law index a is similar for low, intermediate 
and high-mass protostars. In all three cases the average a is 1.5. 
We find that about 60 % of the protostars are well modelled by 
envelopes with 1.5 < a < 2.0. This indicates that 60 % of the 
protostar envelopes are consistent with the standard model of 
free-fall collapse from an initially singular isothermal sphere, 
the so-called inside-out model (Shu 1977). However, it appears 
that 35% of the sample are reproduced by envelopes with 
much smaller indexes, namely 0.5 < a < 1.5. The theoretical 
interpretation of theses low a values is not straightforward. 
The phenomenon was alread y noted in prev i ous studies of 
low-mass protostars ( e.g. lAndre et aLl 119931: IChandler et al.l 
119981: iMotte & Andrei l200ll) . Various hypotheses to explain 
theses low a values have been evoked in the literature. One 
possibility is that the envelopes with low a are described by the 
collapse of an initi ally logptropic sphere, rather than a singular 
isothermal s phere dLizano & Shulll989l: [McLaug hlin & Pudritzl 
1 1 9961 119971: TAndre et al.ll2000b . Basically, the logotropic model 
assumes that the gas pressure across the initial condensation 
depends logarithmically on the density, giving rise to a flatter 
density profile in the static part and at the infall/static interface 
of the envelope. Because the structure of prestellar cores is well 
described by a flat density profile in the inner reg ion (a ~ 0) 
and a power law index of ~ 2 in the outer part (JVisser et al.l 
2002; Andre et al. 2000), it has been mentioned in the past that 
a lower value of a could probe a younger protostar. However, 
later systematic studies have not supported this interpretation 
(e.g. |j0rgensene t al. 2002). Another possible explanation is that 
the envelope is flattened for example because of the presence of 
a magnetic field (e.g. Li & Shu 1996; Hennebelle & Fromang 
2008). However, testing this hypothesis is not trivial as it re- 
quires to solve the radiative transfer in a 2D geometry. A simple 
toy model which assumes constant temperature and optically 
thin emission suggests that a cannot be lower than 1 even in 
the extreme case of a flattened structure with axis ratio of 1:10 
seen face-on. The few cases with a > 2 are easier to explain 
and may be due to t he presence of one or more high-density 
structures, like discs ( J0rgensen et al. 2007), embedded within 
the envelope. 

Envelope mass M env : 

It is very clear from Fig. [T] that more luminous sources have 
larger envelope mass. Hatch ell& van der Takl ([2003) report the 
mass of the envelopes within 1 pc (em pty triangles in Fig. \Tj 
Therefore, the envelope masses from 



Hatchell & van der Tak 
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(2003) are considered as lower limits for M env , and are not taken 
into account in the correlation coefficient computation. The 
correlation coefficient is 0.96 with a probability for a chance 
correlation of ~ 10" 20 , pointing out a strong relation between the 
two variables. Assuming that the luminosity is entirely due to 
the gravitational energy released, the luminosity-envelope mass 
relation suggests that a similar relation exists between the mass 
accretion rate and the mass of the envelope. However, while the 
result that more luminous sources have larger envelope mass 
agrees with theoretical expectations, it should be considered 
with caution. Indeed, the observed correlation between two 
quantities does not necessarily imply a physical correlation. The 
envelope-mass-derivation differs between studies. Furthermore, 
uncertainties in the distance to sources, particularly in the 
high-mass cases (uncertainties of several kpc), and resolution 
limits of distant objects can introduce errors and significant 
observational bias. 

Envelope radius r out : 

Similar to M env , the outer radius of the envelope r out increases 
with increasing luminosity, varying from ~ 10 4 AU for low- 
mass protostars to ~ 10 6 AU for high-mass protostars. The 
correlation coefficient between r out and the luminosity is very 
high, 0.90, with a probability for a chance correlation of 10" 22 . 
Note that we checked that the relation is real and not simply due 
to the distance to the source combined to the angular resolution 
limits. This point is illustrated by the plot of the measured 
angular size, #, of the sources as function of the luminosity in 
Fig. Q] (panel f). The figure shows that 6 is not decreasing with 
luminosity, which excludes the possibility of an observational 
bias due to the limited angular resolution of the observations. 

Average density <n>: 

The average density < n > in the envelope is derived from 
r out and M env for each source. Unfortunately, < n > could be 
derived only for intermediate and some high-mass sources, 
giving an average of about 4xl0 5 and 3xl0 6 cm" 3 , respectively, 
over a spread in luminosity of about four orders of magnitude. 
Although the average < n > is one order of magnitude higher 
in high-mass protostars versus intermediate mass protostars, 
the correlation coefficient between < n > and the luminosity is 
only 0.55. The lower limits for low-mass and the highest mass 
protostars do not allow any firm conclusion, but it seems that 
there is little difference in the average density of these envelopes. 

Density at 1000 AU: 

There is an apparent increase of the density at 1000 AU with 
increasing luminosity, going from about ~ 10 6 to ~ 10 8 cm -3 
for source luminosity varying from 1 to ~ 10 5 L©. There are, 
however, some exceptions on the high-mass side. Indeed, the 
increase of the density at a given distance is consistent with 
the finding of an increasing envelope radius and approximately 
constant average density for the envelope (see the two items 
above). 

Density n\ox at 10 K: 

For a smaller sample, formed by low- and intermediate mass 
protostars only, it is possible to compare the density hiok at 10 
K, which is an indication of the density of the parental cloud 
(the gas shielded by UV photons). Figure [2] shows that n\0K 
is between 2 x 10 4 and 6 x 10 5 cm -3 . This quantity, however, 
is not correlated with the source luminosity (varying by three 
orders of magnitude). Therefore there is no evidence that the 
outer density plays a large role in the determination of the final 



star mass, a rather important and surprising result, which needs 
further confirmation. 

Summary: 

To summarise, the major result of this section is that the proto- 
star luminosity (namely the mass of the final star) seems prefer- 
entially linked to the size (or mass) of the envelope, rather than 
to the parental cloud density, and that most of the envelope ends 
up having a centrally condensed, free-fall density distribution. 
Furthermore, and maybe even more important, there is a con- 
tinuity in the parameters of the envelopes, going from low- to 
high-mass protostars. It appears that there is no important differ- 
ence to the trigger or process of star formation for these two mass 
regimes. The intermediate mass protostars have allowed a bridge 
between the low- and high-mass sources, with no apparent ob- 
servational discontinuity. However, one has to keep in mind that 
these results are based on single dish observations and are there- 
fore driven mainly by the outer region of the envelopes (> 10"). 
A more accurate analysis of the smaller scale structure ( cavities, 
density power law index changes...) will require interferometric 
observations. 



4.2. The problem of the underestimated 24 ^m flux 

As mentioned in the previous sections, our modelling fails to 
reproduce the observed flux at 24/mi by several orders of mag- 
nitude. This is certainly not due to a numerical problem in the 
computation but is rather a real problem: our model misses some 
key element. A possibility evoked in the literature i s the presence 
of a l arge spherical cavity within the envelope ( J0rge nsen et al] 
2005) which could significantly reduce the optical thickness at 
24/mi. 

To check this hypothesis, we carried out a few tests using 
Cep E-mm as a representative case. The results of the tests are 
shown in Fig. [3] First, we added a 1800 AU radius cavity to the 
best- fit model of Table HI As expected, the emission at 24 /mi is 
increased by several orders of magnitude. However, this model 
badly fails to reproduce the brightness profiles (as shown in the 
figure). When a correct procedure is carried out that takes into 
account the variation in the brightness profiles, namely a min- 
imisation of the x^ varying all the envelope parameters, the sit- 
uation returns to the original (best- fit) underestimation of the 24 
/mi flux. In fact, the introduction of a cavity leads to an increase 
of the density power law index, a, to compensate for the result- 
ing flattening of the model brightness profiles. As a result, the 
mass found in the outer envelope is lower so that the the overall 
dust optical depth must be increased to reproduce the integrated 
sub-millimetre fluxes. The dust opacity at 24 /mi again becomes 
large and very little emission is able to escape at these wave- 
lengths. These tests suggest that to solve the 24 /mi flux problem 
it is necessary to have a low-opacity escape route for the 24 /mi 
photons and a thick enough envelope to fit the brightness profiles 
: a simple cavity does no t suffice. 

lLarsson et all (I2Q0QL l2002h similarly found that the 24 //m 
flux is underestimated toward Serpens-FIRS 1 , using different 
tools and observations. In their first study, Larsson et al. used a 
ID radiative transfer code assuming a spherical envelope model 
with a single power law density to reproduce the SED (similar 
to our approach). Their model underestimates the observed flux 
near 24 /mi by several orders of magnitude. In their second study, 
they modelled the envelope with a 2D radiative transfer code, in- 
cluding a biconical cavity. This model reproduces fairly well the 
SED from the mid-IR to millimetre wavelengths, but again un- 
derestimates the observed flux at 24 /mi by about a factor 3-10. 
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Fig. 1. Plot of physical quantities within the envelopes surrounding low-, intermediate, and high-mass protostars as a functio n 
of their bolometric luminosity. Each symbol represents a differe nt group of sources as foll ows: crosses: J0rgens eii et al.1 ([2002); 
filled squares: IShirlev et al.l (l2002l): st ars: this study and| Crimier et al. (2009a.b); diam onds : [W illiams et al. (2005); filled diamonds: 
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represent upper/lower limits. The solid line represents the linear fit of the points obtained when the correlation is significant (Table 
[6]). The dashed line and the dotted line on panel f) represent the angular resolution of the SCUBA maps at 850 and 450 /zm, 
respectively, which were used in all the in the studies reported in this panel. 



The biconical cavity treatment alleviates but does not resolve 
the missing flu x problem. Systemat ic stu dies of high-mass pro- 
tostars by IVan der Tak et all d2000h and IwTlliarns et al] (120051) 
have led to similar results. Willi ams et al.l (2005) modelled 36 
high-mass protostellar objects at 850 //m, using DUSTY, as in 
the present study. The majority of their best-fits fail to repro- 
duce the flux around 24 /mi. They discuss the possibility of the 
contribution of high accretion rates, which would significantly 
increase the near-IR flux density (Osorio et al. 1999). Another 
contribution could come from the presence of circumstellar disks 
or the stochastic heating of small grains, which would alter the 
emission of the envelope and produce more short wavelength 
photons (e.g. lSellgren et alJll983HDraine & Lill200ll) . However, 
the envelope is optically thick at these wavelengths, making the 
emission very sensitive to deviations from the assumed spherical 
shape. These additional processes still require a low-opacity es- 
cape route to exists (also mentioned by 1 Van der Tak et al]|2000h 
to explain the missing 24 /mi flux, for example a biconical 
cavity excavated by the outflow. Note that all the intermediate 
mass protostars studied here are associated with outflows, except 
OMC2-FIR4 for which t he model indeed fits the observed flux 
at 24//m (see ICrimier et al. 2009a). Finally, Van Der Tak et al. 
also suggest the possible evaporation of grain ice mantles close 



to th e star, which would decrease the 20 //m optical depth by 
30% (lOssenkopf & Henning|[T994l) in the T dust > 100 K region. 



5. Conclusions 

We have derived the physical structure of the envelopes of five 
IM protostars, with luminosities between 30 to 1000 L©. The en- 
velope dust density and temperature profiles were determined by 
means of the ID radiative transfer code DUSTY, using all con- 
tinuum observations from the literature. The analysis assumed 
that the density profiles follow a single index power law and ob- 
tained self-consistently the temperature profile. The best-fit en- 
velope models well reproduce the observations, namely the sub- 
millimetre radial brightness profiles and the SED between ~60 
yum and 1.3 mm, for each source. However, the model underes- 
timates the 24//m emission by several orders of magnitude. We 
ran test models to better understand what ingredient is missing 
and conclude that a "simple" cavity is not enough to reproduce 
the 24/mi observations. Apparently, the missing ingredient is a 
low-opacity escape route plus a warm dust contribution inside 
the envelope (circumstellar disc, warm outflow-excavated cav- 
ity...). 

The gas density and temperature profiles were derived by as- 
suming a constant dust-to-gas ratio and by computing the gas 
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Fig. 2. Plot of the density at 10 K as function of the bolometric luminosity of a sample of low- and intermediate mass protostars. 
Data are from J0rgensen et al. (2002), (crosses), and the present study, (stars). 



thermal balance at each point within the envelope. Because the 
gas equilibrium temperature strongly depends on the water abun- 
dance in the interiors of the envelopes, we also computed the ex- 
pected water emission for each source. We found that the gas and 
dust are thermally coupled across the envelope with differences 
less than 5 K in three out of five sources. In IC1396 N BIMA 2 
and NGC7129 FIRS 2 the gas is colder than the dust by at most 
40 K, in a small region just where the icy mantles are predicted 
to sublimate. The predicted water line fluxes are consistent with 
the upper limits derived by the ISO observations. 

One of the major goals of the present study was to "use" the 
IM protostars as a bridge between the low- and high-mass pro- 
tostars with the hope that this will aid our understanding of the 
star formation process at either end. When comparing the char- 
acteristics derived by the modelling of the envelopes of low-, 
intermediate, and high-mass protostars, it appears that there is a 
smooth transition between the various groups. This suggests that 
there are basically no different triggers or processes between 
these mass regimes. The power law index a is similar in all three 
groups of objects. The majority of the sources have a between 
1.5 and 2. This is consistent with the theory of isothermal col- 
lapse from an initially singular isothermal sphere, the so-called 
inside-out expansion- wave collapse (Shu 1977). Regardless of 
the mass group, a few sources have a lower than 1.5, point- 
ing perhaps to the collapse of an initially logotropic, virialised 
sphere (Lizano & Shu 1989, McLaughlin 96, and McLaughlin 
97). Finally, the luminosity (mass) of the star depends on the 
size of the envelope, but does not depend on the density at a 
given temperature (for example at 10 K or 100 K). 
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Fig. 3. Cep E-mm physical structure, brightness profiles, and SED from the best-fit model (black solid lines) and after adding an 
1800 AU radius cavity to the best-fit model (red dashed-dotted lines). The dust temperature and H2 density profiles are reported 
in the upper left panel and upper right panel, respectively. The observed brightness profiles at 450 yum, 850 yum and the SED, 
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best-fit model. 
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Appendix A: Results for individual sources 

A.1. CB3-mm 

AAA. Source background 

The CB3 Bok globule is located at ~ 2.5 kpc 
(lLaunhardt & Henningl Il997l: IWang et al] 1 1995b on the near 
side of the Pe r seus arm of the Galaxy. Using ISO data, 
lLaunhardt et all (1 19981) carried out a multi-wavelength study 
of the CB3 globule and derived a total mass and bolometric 
luminosity of the order of 400 M© and ~ 1000 L©, respec- 
tively, for the entire globule. The globule hosts ~ 40 NIR 
sources, 22 of which are likely low- and intermedi ate mass 
proto st ars in different stages of evolution (Launh ardt et al] 
119981: lYun & Clemenl I1995L 11994]) . CB3-mm, the bright- 
est millimetre source in the globule, was first detected by 
lLaunhardt & Henning! (1 1997b a nd subsequently observed in 
the sub-millimetre by Hua rd et all (l2000h. The hig h luminosity 
of the source evaluated by Launhardt & Henning (1997), L 
~ 900 L©, suggests that CB3-mm is a n intermediat e mass 
object. A recent interferometric study by iFuente et all (2007) 
showed that CB3-mm contains t wo compact cores at 3 mm 
separated by 0.3 pc (~ 0.43"). lYun & Clemensl (11994 also 
detected a molecular bipolar outflow in CO, elongated in the 
northeast- southwest direction, as sociated with H2O masers 
(Ide Gregorio-Monsalvo et al] [2006). This outflow has been 
mapped in various molecular lines by Codella & Bachiller 
(1999), who concluded that it originates from CB3-mm. The 
same authors concluded that CB3-mm is probably a Class 
source. In the present study we re-evaluated the luminosity of 
this object, L ~ 1000 L©, making CB3-mm an intermediate 
mass protostar. 

A.1.2. Analysis 

The continuum maps used for the CB3-mm analysis are pre- 
sented in Fig. IA.1I The maps at 450 /mi and 850 fim were ob- 
tained on 1997 December 18 as a part of project m96bi28 and on 
1998 August 10 as a part of project m98bc21, respectively. The 
flux profiles obtained at each wavelength are shown in Fig. IA.2l 
The Spitzer observations were obtained on the 20th September 
2004 as part of the programme "Comparative Study of Galactic 
and Extragalactic HII Regions" (AOR: 63, PI: James R. Houck). 
The integrated fluxes used for the analysis are reported in Table 
[2]and in Fig. IA.21 The fluxes at 70 yum, 450 yum and 850 /mi were 
obtained by integration over a ~ 50"radius. The uncertainty el- 
lipse position of the IRAS observations are reported on the 450 
yum maps in Fig. IA.U Note that the map at 24 /mi presents two 
sources separated by - 12", centred on the sub-millimetre source 
emission (see Fig. IA.ll) . Therefore, the flux 1.4 ± 0.8 Jy at 24 yum 
was obtained by adding the integrated flux over the two sources 
separately. We also report on the SED the lower limits from 
Plateau de Bure (PdB) interferometric fluxes at 1.3 mm and 3 
mm (Fue nte et al] |2007). 

A. 1.3. Best-fit 

Table [4] presents the set of parameters a, Y, and T100, which best 
reproduce the observations, and sum mari ses some relevant phys- 
ical quantities of the model. Figure \AJ2\ shows the best derived 
brightness profiles and SED against the observed ones, while 
Fig. IA.3I shows the^ 2 contour plots. The dust density and tem- 
perature profiles of the best-fit model are reported in Fig. IA.4I 
Although the observed flux profiles and SED fluxes from 60 yum 
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Fig. A.l. Continuum emission maps around CB3-mm at 24 /mi (upper left panel), 70 yum (upper right panel), 450 /mi (lower left 
panel), and 850 yum (lower right panel). The grey scale contours mark the continuum flux from 10 % to 90 % of the peak emission 
by steps of 20 %. The uncertainty ellipse position of the IRAS observations toward CB3-mm are represented by the ellipse in the 
450 jim map. 



to 850 yum are well reproduced by the model, the observed flux at 
24 //m is underestimated by about two orders of magnitude. The 
model agrees well with the lower limits fr om Plateau de Bur e 
interferometric fluxes at 1.3 mm and 3 mm dFuente et al .1120071) . 
however. 

The best-fit model, obtained with an a of 2.2, is the steep- 
est among the five sources modelled in the paper and leads to a 
strong gradient in the density profile, from ~ 10 9 to ~ 10 3 cm -3 
with a density and radius at 100 K of 7.5xl0 7 cm -3 and 700 AU, 
respectively. 

Although CB3-mm is the largest and brightest of the five 
sources investigated in this paper, its distance of 2500 pc makes 
it the least resolved source. Ind eed, t he brightness profiles at 
each wavelength presented in Fig. lA.2l are close to the beam pat- 
terns of the telescope, showing that the source is barely resolved. 
Consequently, the value of the envelope radius r out suffers a rel- 
atively large uncertainty. Considering \h&x 2 re d contour at 10% of 
the minimum x 2 ed (see Fig. IA.3b in order to estimate the uncer- 
tainty, we obtain r out = (1.0 ±3' 3 )xl0 5 AU. On the contrary, the 
envelope column density is quite insensitive to r out , due to the 
high power law index of the density profile. Finally, in minimis- 
ing the xl ed > we a ^ so var i e d the source luminosity from 800 L© 
to 1200 L© and found that the best-fit is obtained with source 
luminosity equal to 1000 L . 



A.2. Cep E-mm 

A.2.1 . Source background 

Located in the Cepheus OB3 association at a distance of 730 
pc (Sargent 1977; Crawford & Barnes 1970), Cepheus E is the 
second most massive and dense clump of this region (Few & 
Cohen 1983). Cepheus E hosts the source IRAS 23011+6126, 
which is assoc iated with Cep E- mm, catalogued as a Class 
protostar by iLefloch et al.l {1996). Targeted by several contin- 
uum andfines_studies, Cep E-m m has been observed with IRAS 
(Palla et alJ fl993b. IRAM 30m (ILefloch et al II1996I: IChini et al.l 



200 lb. SCU BA (IChini et al.ll2Q0lh ISO (Froebrich et al.ll2003l). 
and Spitzer (Noriega-Crespo et al. 2005). All these studies con- 
firm the Class status of Cep E-mm and constrain the total mass 
and bolometric luminosity of the source to 7-25 Mq and ~ 80- 
120 Lq, respec tively. Finally, a bipolar molecular outflow, first 
reported by iFukui et all (Il989l) . is associated with Cep E-mm. 
The pro perties of the outflow have been thoroughly analysed 
by Eislof fel etal.1 (119961) lAyaja et al.l (l2000h . lMoro-Martin et al.l 
(20 0lh. and ISmith et al.l i (l20Q3h. The H 2 and [Fell] study by 
Eisloffel etal. I (Il996l) shows a quadrapolar outflow morphol- 
ogy, suggesting that the driving source is a binary. This outflow 
morpholog y has been confirmed by sub-mm and near-IR ob- 
servations (Ladd & Hodaprj fl997), and by Spitzer observations 
(iNoriega-Crespo et al.ll2005l) . 
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Fig. K3.x 2 re d contour pl° ts (Y,#) for the CB3-mm model. In these computations, Tioo is equal to 5.8. The contours show the loci 
oi\hQX 2 re d vames equal to 1.1, 1.2, 1.5 and 2 times the minimum x 2 red (10, 20, 50 and 100%). The cross shows the position of the 
best-fit. 



A.2.2. Analysis 

The continuu m ma ps used for the Cep E-mm analysis are pre- 
sented in Fig. IA.5I The maps at 450 /mi and 850 jum were ob- 
tained in August 1997 as a part of the project m97b u87 and the 
observations are described in detail by IChini et al.l (1200 ll) . The 
flux profiles obtained at each wavelength are shown in Fig. IA.6l 
The Spitzer observations were obtained on the 2003 September 
29 as part of the programme "MIPS/IRAC imaging of protostel- 
lar jet HH 212" (AOR: 1 063, PI: Alberto Noriega-Crespo) an d 
are described in detail by Norieg a-Crespo & Garnavich (2001). 
The integrated fluxes at 450 /mi, 850 urn, an d 1300 /mi were 
retrieved from the literature (IChini et al.ll200ll) . The uncertainty 
ellipse position of the IRAS observation is reported on the 450 



yum map in Fig. IA.5I The fluxes at 24 fim and 70 yum were ob- 
tained by integration over a ~ 30"and ~ 60"radius, respectively. 
Note that to compute the flux uncertainties, we varied the inte- 
gration radius of about 50% and found variations of < 30%. The 
integr ated fl uxes used for the analysis are reported in Tableland 
in Fig. El 



A.2.3. Best-fit 

Table |4] presents the set of parameters a, Y, and tioo, which best 
reproduce the observations, and su mma rises some relevant phys- 
ical quantities of the model. Figure lA!6l shows the derived bright- 
ness profiles and SED against the observations. The dust density 
and temperature profiles of the best-fit model are reported in Fig. 
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Fig. A.4. CB3-mm dust temperature (left panel) and H2 density (right panel) profiles from the best-fit model. The dashed line 
represents the difference between the dust and gas temperature derived using X(H20); n = lxlO -5 and X(H20) 0Ut = lxlO -7 (scale 
on right axis). 
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Fig. A.5. Continuum emission maps around Cep E-mm at 24 yum (upper left panel), 70 /mi (upper right panel), 450 yum (lower left 
panel), and 850 yum (lower right panel). The grey scale contours mark the continuum flux from 10 % to 90 % of the peak emission 
in steps of 20 %. The uncertainty ellipse position of the IRAS observation toward Cep E-mm is reported on the 450 /mi map. 



IA.7I Although the observed flux profiles and SED fluxes from 60 
yum to 1300 yum are well reproduced by the model, the observed 
flux at 24 yum is underestimated by about ~3 orders of magni- 
tude. Similarly to CB3-mm, the envelope model derived for Cep 
E-mm with an a of 1.9 is very peaked, leading to a strong gradi- 
ent in the density profile, from ~ 10 9 to ~ 10 4 cm -3 with a den- 
sity and radius at 100 K of 2.0xl0 8 cm -3 and 223 AU, respec- 
tively. Considering the^^ contour at 10% of the minimum x 2 red 
in order to estimate the uncertainty, we obtained r out = (3.6+0.7) 
xlO 4 AU. Finally, in minimising the^ £D , we also varied the 



source luminosity from 70 L© to 130 L© and found the best-fit 
for a source luminosity equal to 100 L©. 

A.3. IC1396NBIMA2 
A.3.1 . Source background 



IC1396 N is a bright globule located 750 pc (lMatthewslll979b 
from t he Sun, near the border of the IC139 6 extended HII 
region (lOsterbrockl Il957l: IWeikard et aD Il996l) and at a pro- 
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jected distance of 11 pc north of the 06.5 star HD 206267 
which ionises the region. IC1396 N is associated with the 
source IRAS21391+5802. I ts strong submillimetre and mil- 
limetre continuum emission (JWilking et al.lll993t [Susitani et al. | 
12000b ICodella et aD l200ll) . high-density gas (Serabvn etalJ 



1993 



1995 



ICesaronietaD 119991: [Codella et al. 2001; Beltra n et al 



2004) , and water mase r emission (JFelli et al.lll992UTofani et al. 






[ iPateletaLlbOOQl) reveal that IC1396 N is an active site of 



star formation. The bolometric luminosity is estimated to range 
from 235 L Q dSaraceno et al.l H996) to 440 L Q ( Sugitaniet all 
I2000h . 

Using B I MA interferometric millimetre observations, 
iBeltran et all (120021) detected three sources (BIMA 1, BIMA 
2 and BIMA 3) deeply embedded in the globule. BIMA 2 has 
the strongest millimetre emission and is also the most massive 



object. The authors concluded that BIMA 2 is most likely an IM 
protostar, while BIMA 1, and BIMA 3, located at ~15"west and 
south-east respective ly, are less massi ve a nd/or more evolved . 
The recent studies bv lNeri et all (120071) and lFuente et al.1 (l2007h . 
using Plateau de Bure interferometric observations, provide a 
detailed study of BIMA 2 and showed that BIMA 2 contains at 
least three dense cores separated by ~1". 

Finally, the region presents an extended (~3 arcmin) CO 
bi polar outfl ow (Sugitan i et al.l 1 1989b which has been ma pped 
by ICodella et alJ (l200lh . The study bv iNisini et aD (l2001h . re- 
vealed several H2 jets inside the region. A more recent analysis 
(IBeltran et alJ l2002. 2004) shows that BIMA 1 and BIMA 2 are 
associated with north-south and east- west bipolar molecular out- 
flows, respectively. 
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Fig. A.8. Continuum emission maps around IC1396 N BIMA 2 at 24 /mi (upper left panel), 70 /mi (upper right panel), 450 /mi 
(lower left panel), and 850 /mi (lower right panel). The grey scale contours mark the continuum flux from 10 % to 90 % of the peak 
emission by steps of 20 %. The positions of BIMA 1, BIMA 2 and BIMA 3, as well as the uncertainty ellipse position of the IRAS 
observations are reported on the 24 /mi map. The dashed regions on the 450 //m and 850 /mi maps mark the regions excluded to 
derive the flux profiles. 



A.3.2. Analysis 

The continuum maps used for the IC1396 N BIMA 2 analysis are 
presented in Fig. IA.8l The SCUBA maps were obtained in March 
2002 as a part of the project m02au07. The brightness profiles 
(see Fig. IA.9I) are derived excluding the regions contaminated 
by the presence of BIMA 1 and BIMA 3 (dashed regions in Fig. 
IA.8D . The Spitzer observations were obtained in October 2004 
as part of the programme "Star Formation in Bright Rimmed 
Clouds" (AOR: 202, PI: Giovanni Fazio). The integrated fluxes 
used for the analysis are reported in Table [2] and in Fig. IA.9I 
The integrated fluxes at 450 /mi and 850 //m were obtained by 
integration over a ~ 35". We varied the integration radius from 
20"to 90", and found variations on the fluxes < 30%. The large 
uncertainty in the integrated flux at 24 /mi is due to the proximity 
of BIMA 3 (see 24 /mi map in Fig. IA.8I) . Note that due to the 
large IRAS beams at 60 //m and 100 /mi, BIMA 3 introduces 
significant spurious signal. We thus used the IRAS and 70 /mi 
Spitzer fluxes as upper limits on the SED (Fig |A.9l) . 



A.3.3. Best-fit 

Table |4] presents the set of parameters a, Y, and rioo, which 
best reproduce the observations, and su mmar ises some relevant 
physical quantities of the model. Figure [A79l shows the relevant 
derived brightness profiles and SED against the observations. 



The dust density and temperature profiles of the best-fit model 
are reported in Fig. IA.10I The observed flux profiles and SED 
fluxes at 450 /mi and 850 /mi are well reproduced by the model. 
Considering the IRAS and 70 /mi Spitzer fluxes as upper limits 
(see Sect. IA.3.21) , the minimum Xsed * s obtained for a source 
luminosity equal to 150 L©, similar to the value suggested by 
iBeltran et al.l (120021) . The observed flux at 24 //m is underesti- 
mated by about one order of magnitude. Considering the x^ ed 
contour at 10% of the minimum x 2 red to estimate the uncertainty, 
we obtained r out = (3.0+0.6) xlO 4 AU. Note that we obtain sim- 
ilar results without subtracting the dashed region of Fig. IA.81 
confirming the weak contribution of BIMA 1 and 3 to the sub- 
millimetre emission. 

Finally, contrary to the previous sources, the gas tempera- 
ture profile obtained for IC1396 N BIMA 2 is slightly decou- 
pled from the dust temperature in the inner region (see Fig. 
lA.lOb . The origin of this thermal decoupling is discussed in de- 
tail in[Crimier et al. (2009a). Briefly, the decoupling is due to a 
competition between the water abundance (the main coolant of 
the gas in the inner region) and the dust density, which mainly 
regulates the balance energy in the inner part of the envelope. 
This thermal decoupling is stronger in IC1396 N BIMA 2 (and 
NGC7129 FIRS 2) than in the other sources because of the lower 
density in the inner region. The thermal decoupling reaches a 
maximum value of about 40 K (-15% of T^ust) for an abun- 
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dance of X(H20)/„ = 1x10 5 and about 10 K for an abundance 
ofX(H 2 O) m = lxl0" 6 . 

A.4. NGC7129 FIRS 2 

AAA. Source background 

NGC7129 is a reflecti on nebula located in a complex and 
active molecular cloud (lHartigan & Ladal Il985l; iMiranda et al.l 
Il993l) and estimated to be at a distance of 1250+50 pc 
(Sh evchenko & Yakubov|[l989l) from the Sun. The region con- 
tains several Herbig AeBe stars, which illuminate the nebula. 
NGC7129 FIRS 2 has been detected in the far-infrared by 
iBechis et al.l dl978h and lHarveyetaD (11984 . FIRS 2 is not de- 



tected at optical or near-infrared wavele ngths. Its position coin- 
cides with a 13 CO column density peak (iBechis et al.lll97 8) and 
a high-density NH3 cloudlet (Guesten & Marcaide 1986). It is 
also close to an H 2 maser (Rodriguez et al.lll980h . NGC7129 
FIRS 2 has been classified as an IM Class source by lEiroa et al.l 
(1998), who carried out a multi- wavelength study of the contin- 
uum emission from 25 yum to 2000 jim. These authors estimate 
a total mass and bolometric luminosity of ~ 6 M© and ~ 430 
L©, respectively. In addition, Edwards & Snell (1983) detected 
a bipolar C O outflow associated with FIRS 2. The interferomet- 
ric study of iFuente et al.l (1200 ll) has confirmed this bipolar out- 
flow and pointed out a quadrapolar morphology to the flow. This 
quadrapolar morphology seems to be due to the superposition 
of two flows, FIRS 2-out 1 and FIRS 2-out 2, likely associated 
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with FI RS 2 and a more evo lved star (FIRS 2-IR), respectively. 
Finally. iFuente et al. (2005 a. b) carried out an extensive chemical 
study of FIRS 2 providing the first detection of a hot core in an 
IM Class 0. Based on all these studies, FIRS 2 is considered the 
youngest IM object known at present. 



A.4.2. Analysis 

All the data and continuum maps used for the NGC7129 FIRS 
2 analysis are reported in Table [2] and Fig. IA.1U T he maps 
at 450 yum and 850 yum were obtained on the 1998 August 26 
as a part of the project M98BU24. We also used the flux pro- 
file and the integrated flux at 1300 //m extract ed from t he sin - 
gle dish map at 1300 jim, presented in Fue nte et all (120011) . 
Following iGreve et all (Il998l) . the beam at 1300 /zm is assumed 
to be a combination of three Gaussian curves, with HPBWs of 
11", 125", and 180", with amplitude ratios of 0.975, 0.005, and 
0.001, respectively. The Spitzer observations were obtained on 
the 2003 September 19 as part of the programme "Protostars 
and Proto-Brown Dwarfs in a Nearby Dark Cloud" (AOR: 34, 
PI: Tom Megeath). The integrated fluxes reported in Tableland 
in Fig. IA.12I were obtained by integration over a ~ 40"radius. 
The IRAS integ rated fluxes at 60 yum and 100 yum were taken 
from lEiroa et aD (119981) . 



A.4.3. Best-fit 

Table |4] presents the set of parameters a, Y, and rioo, which best 
reproduce the observations, and su mmaris es some relevant phys- 
ical quantities of the model. Figure IA7T21 shows the relevant de- 
rived brightness profiles and SED against the observations. The 
dust density an d tem perature profiles of the best-fit model are 
reported in Fig. IA.13l The observed flux profiles and SED fluxes 
from 60 yum to 850 yum are well reproduced by the model. In 
minimising the^ £D , we also varied the source luminosity from 
400 L to 600 L© and found the best-fit for a source luminosity 
equal to 500 L . The observed flux at 24 yum is underestimated 
by about two orders of magnitude. Considering the^^ contour 
at 10% of the minimum x 2 red to estimate the uncertainty, we ob- 
tained r out = (1.9+0.2) xlO 4 AU. 

Finally, similarly to IC1396 N BIMA 2, the gas temperature 
profile obtained in NGC7129 FIRS 2 is slightl y deco upled from 
the dust temperature in the inner part (see Fig. IA.13b . 



A.5. Serpens FIRS 1 

A.5.1. Source background 

Since its recognit ion as a very active star forming region 
(iStromet al.lll974l) . the Serpens molecular cloud has been the 
target of many observational studies. Several studies have been 
dedicated to evaluate the distance of this cloud from the Sun, 
yielding a rang e of distan ces from ~ 210 to 440 pc. Following 
the review by lEiroa et al.l (120081) . which reports and compares 
the distance results from all the different studies, we adopt 
a distance of ~ 230+20 pc. Serpens is a young protoclus- 
ter whose members are in many different evolutionary stages. 
About 100 embedded YSQ, protos t ars, and prestellar clumps 
(Strom e taD 119761: lEi roa & Casali 119921; ICasali et al.l 119931: 
Hurt&Barsonv 1996; Giovannetti et al. 1998; Testi & Sargent 
1998; Kaas 1999) , as well as molecular outflow s ( White etaD 
1995; Wolf-Chas e"etaDll998l:lwriliams & Mversll2000h and hT 



emission (e.g. lEiroa et al.l 119971: iHerbst et all 1 19971: iHodappI 
1999) have been detected in this region. 

Serpens FIRS 1 is located near the centre of the Serpens main 
core and is the most luminous object embedded in the cloud. 
Several continuum studies classify it as Class source with a 
bolometric luminosit y estimated to range from 46 Lp to 84 Lq 
(Ha rvev et all 1 19841: ICasali et al] 1 19931: iHurt & Barsonvl 1 19961: 
lLarsson et alj|2000l) . The latest estimation of its luminosity sug- 
gests that Serpens FIRS 1 is on the low/intermediate mass pro- 
tostar border. The envelope physical structure of Serpens FIRS 1 
has been modelled by Hogerheijde et al. (1999) using millime- 
tre interferometric observations of the continuum and molecular 
line observations. The Serpens FIRS 1 SED was first modelled 
by means of a sph erical envelope I D radiative transfer c ode by 
Lar sson et al.l (2000). Subsequently, Larsson et al.l (2002) used a 
2D radiative transfer code and a torus model to reproduce the 
SED and line observations from ISO. FIRS1 drives a molec- 
ular outflow, which is orientated at a position angle of 50 de- 
grees, coincident with a triple radio source with symmetric lobes 
dRodriguezetaDl989h . 

A.5.2. Analysis 

The continuum maps used for the Serpens FIRS 1 analysis are 
presented in Fig. IA.14I The SCUBA maps were obtained in 
January 1998 as a p art of the project m97bc30. The flux pro- 
files (see Fig. IA.15I) are derived excluding the dashed regions 
in Fig. lA.14l to avoid the non-spherical extended emission. The 
Spitzer observations were obtained in October 2004 as part of 
the programme "From Molecular Cores to Planets, continued" 
(AOR: 174, PI: Neal Evans). The integrated fluxes used for the 
analysis are reported in Table [2] and in Fig. IA.15I The inte- 
grated fluxes from 70 /im to 850 yum were obtained by integration 
over a ~ 40". We also report on the SED the lower limits from 
Plate au de Bure (PdB) interferometric fluxes at 1.3 mm and 3 
mm dFuente et alJl2007h . 



A.5.3. Best-fit 

Table [4] presents the set of parameters a, Y, and rioo, which 
best reproduce the observations, and sum marise s some relevant 
physical quantities of the model. Figure IA.15I shows the rele- 
vant derived brightness profiles and SED against the observa- 
tions. The dust density a nd tem perature profiles of the best-fit 
model are reported in Fig. IA.16l The model well reproduces the 
observed flux profiles and SED fluxes from 70 /mi to 850 /im 
and agrees well with the lower limits from Plateau de Bure in- 
terferometric fluxes at 1.3 mm and 3 mm dFuente et al.ll2007h . In 
minimising the xIed> we also varied the source luminosity from 
20 L© to 60 L© and found the best-fit for a source luminosity 
equal to 33 L©. Note that the estimated value of the bolomet- 
ric luminosity is much lower t han the values derived in previ- 
ous s t udies toward this source (lHarvev et al .11 19841 : ICasali et all 
119931: IHurt & Barsonvl 1995 lLarsson et alJl2000h . This is due to 
the adopted smaller distance to the source. Considering \h&x 2 re d 
contour at 10% of the minimum x 2 red to estimate the uncertainty, 
we obtained r out = (5. 9 +0.3) xlO 3 AU. 

Hogerhe ijde et al.l (Il999h studied the envelope physical 



structure of Serpens FIRS 1 using interferometric continuum ob- 
servations at OVRC0 at A = 3.4-1 A mm and single-dish line 
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Fig. A.ll. Continuum emission maps around NGC7129 FIRS 2 at 24 /mi (upper left panel), 70 /mi (upper right panel), 450 /mi 
(middle left panel), 850 /mi (middle right panel), and 1300 /mi (lower centred panel). The grey scale contours mark the continuum 
flux from 10 % to 90 % of the peak emission by steps of 20 %. 



observations from JCMT and CSC0. They adopted a bolomet- 
ric luminosity of 77 L© and fixed the outer r adius of the enve- 
lope R out = 8000 AU for a distance of 400 pc. Hog erheijde et al.l 
(1999) derived an envelope with density profile index a = 2+0.5 
in good agreement with our results. By scaling their results to the 
same distance adopted here (230 pc) their adopted luminosity, ~ 
25 L , and the envelope radius R ou t, ~ 4600 AU, are relatively 
consistant with our results (Table |4]) : L ~ 33 L , R out ~ 5900 
AU. Also their derived gas temperature profile slightly differs 
from ours: for example they obtained T gas ~ 40 K at 200 AU 



Caltech Submillimetre Observatory 



while we obtained T gas ~ 65 K. lLarsson et a l. (2000) studied the 
SED of Serpens FIRS 1, derived from the ISO LWS maps and 
photometric data between ~ 1 and 10 5 /mi retrieved from the lit- 
erature. These authors used i) a single temperature/optically thin 
approximation and ii) a radiative transfer model to analyse the 
SED. In their models they adopted a distance of 3 10 pc, a source 
diameter of 7" and a power law density profile index, a =2. They 
succeeded in reproducing the SED between 40 /mi and 3 mm. By 
scaling their results to the same distance adopted here (230 pc), 
they obtain a luminosity of ~ 40 L , against our estimate of 33 
L , and an envelope mass of 1.1-3.3 M , against our estimate of 
5 M . However, as noticed by the authors themselves, their mass 
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Fig. A.12. Observed brightness profiles at 450 /mi (upper left panel), 850 yum (upper right panel), 1300 yum (lower left panel), and 
the SED (lower right panel) toward NGC7129 FIRS 2. The solid lines report the best-fit model (Table [4]). The dashed lines represent 
the beam pattern of the telescope at 450 yum, 850 yum, 1300 /mi. The vertical dotted lines show the outer radius of the best-fit model. 
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Fig. A.13. NGC7129 FIRS 2 dust temperature (left panel) and H2 density (right panel) profiles from the best-fit model. The dashed 
line represents the difference between the dust and gas temperature derived using X(H20)/ n = lxlO -5 and X(H20) 0Mr = lxlO -7 
(scale on right axis). 

estimate is rather a lower limit, given the optically thin approxi- 
mation they used. 

Finally, the observed flux at 24 /mi is underestimated by 
about two o rders of magn i tude. This phenomenon is also ob- 
served in the lLarsson et al.l (I2002i l2000) studies. We discuss the 
situation at 24 yum in more detail in the the following section, 
Sect. [4] 
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Fig. A.14. Continuum emission maps around Serpens FIRS 1 at 24 yum (upper left panel), 70 jim (upper right panel), 450 yum (lower 
left panel), and 850 yum (lower right panel). The grey scale contours mark the continuum flux from 10 % to 90 % of the peak 
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profiles. 



Source 


CB3 


CepE 


IC1396N 


NGC7129 


Serpens 




mm 


mm 


BIMA2 


FIRS 2 


FIRS1 


RA (2000) 


00:28:42.1 


23:03:12.7 


21:40:41.8 


21:43:01.5 


18:29:49.8 


Dec (2000) 


56:41:59.4 


61:42:27.4 


58:16:13.5 


66:03:25.0 


01:15:18.4 


Distance (pc) 


2500 


730 


750 


1250 


230 


Luminosity (L©) 


1000 


100 


150 


500 


33 



Table 1. The coordinates, the distances and the bolometric luminosities of the sources of the sample. 



Source 


CB3 


CepE 


IC1396 N 


NGC7129 


Serpens 




mm 


mm 


BIMA2 


FIRS 2 


FIRS1 


F 2 4 yU m (Jy) - Spitzer 


1.4 + 0.8 


6.0 + 3.0 


2.3 + 2.0 


6.0 + 3.0 


3.0+1.0 


F 6 (w (Jy) - IRAS 


31 + 15 


61+30 


- 


- 


- 


Fvo^m (Jy) - Spitzer 


58 + 20 


100 + 40 


- 


110 + 60 


135 + 50 


Flo**,, (Jy) - IRAS 


110 + 55 


110 + 55 


- 


- 


- 


F 45 o^ (Jy) - JCMT 


30+10 


44 + 10 


100 + 30 


80 + 40 


130 + 40 


F 85 (w (Jy) - JCMT 


2.0 + 1.0 


4.1 + 1.3 


17 + 7.0 


4.8 + 2.0 


14 + 4.0 


Fisoo^m (Jy) - IRAM30m 


- 


1.0 + 0.2 


- 


2.2 + 0.5 


- 



Table 2. Continuum fluxes integrated over the envelope of each source and used in the SED analysis. 
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Fixed input parameters 



Stellar temperature, T+ 

Dust temperature at r in , T- m 

Dust opacity (0H5) at 100//m, k 10 o 



5000 K 

300 K 

86.5 cmV 1 



Best-fit parameters & physical quantities 



Source 


CB3 


Cep E 


/CI 396 N 


NGC1129 


S erpens 


OMC2 




mm 


mm 


BIMA2 


FIRS 2 


FIRS 1 


FIR4 


Distance, (pc) 


2500 


730 


750 


1250 


230 


450 


Luminosity, (L©) 


1000 


100 


150 


500 


33 


1000 


Dust optical depth at 100//m, rioo 


5.8 


5.0 


1.4 


2.3 


3.0 


0.6 


Density power law index, a 


2.2 


1.9 


1.2 


1.4 


1.5 


0.6 


Envelope thickness, Y (= r out /r in ) 


400 


500 


630 


180 


200 


120 


X^-maps 


0.43 


0.11 


0.55 


0.21 


0.20 


1.2 


^L-SED 


0.85 


0.80 


1.2 


0.87 


4.4 


0.55 


Inner envelope radius, r in , (AU) 


260 


70 


50 


100 


30 


100 


Outer envelope radius, r out , (AU) 


103000 


35800 


29600 


18600 


5900 


12000 


Radius at T dust = 100 K, r 100K , (AU) 


700 


223 


180 


373 


102 


440 


H 2 density at tiook, n , (cm -3 ) 


7.5 x 10 7 


2.0 x 10 8 


4.3 x 10 7 


4.4 x 10 7 


2.2 x 10 8 


4.3 x 10 6 


Envelope mass, M env , (M©) 


120 


35 


90 


50 


5.0 


30 


\*-dust~ *- gasjmaxi v-^J 


6 


5 


40 


40 


10 


150 



Table 4. Summary of the dust radiative transfer analysis for the five IM protostars. The top section lists the fixed input parameters, 
the lower section reports the best-fit parameters, and some relevant physical quantities corresponding t o the best-fit model. The 
last column reports the same parameters for the IM protostar OMC2-FIR4, found with a similar analysis dCrimier et al.ll2QQ9ab . for 
comparison. 



Source 


CB3 


CepE 


IC1396 N 


NGC7129 


Serpens 




mm 


mm 


BIMA2 


FIRS 2 


FIRS1 


A 


Transition 


Model LWS 


Model 


LWS 


Model 


LWS 


Model 


LWS 


Model 


LWS 


Qim) 


3k-k+- 3'k'_k' + 




















ortho 






















75.38 


321 - 2i2 


0.3 


0.2 


1.4+0.3 


0.08 


<3.1 


0.2 


<1.7 


0.5 


2+0.4 


108.07 


22i- lio 


0.1 


0.1 


1.8+0.5 


0.07 


<2.3 


0.1 


<1.5 


0.4 


1.5+0.4 


113.54 


4l4 - 3()3 


0.1 


0.1 


2.2+0.5 


0.04 


<1.2 


0.08 


<1.0 


0.2 


2.9+0.9 


174.62 


3o3- 2i2 


0.06 


0.07 


1.6+1 


0.06 


<2.6 


0.08 


<0.8 


0.2 


2+0.2 


179.53 


2l2~ loi 


0.09 


0.1 


2.9+0.3 


0.2 


<2.2 


0.2 


<1.1 


0.5 


1.4+0.2 


para 






















89.99 


322~ 2n 


0.2 


0.1 


1+0.2 


0.04 


<2.4 


0.1 


<1.5 


0.3 


2.4+0.6 


100.98 


220- In 


0.1 


0.1 


<0.5 


0.04 


<1.0 


0.09 


<0.6 


0.3 


2.6+0.5 



Table 5. The brightest water lines predicted by the our model considering X(H20); W 
with the ISO LWS observations. The upper limits refer to 3<x. All fluxes are in 10" 12 



lxlO- 5 andX(H 2 OW 
erg s" 1 cm" 2 units. 



1x10 7 compared 



Physical parameter 


Linear correlation 


Number of 


Probability for 


Linear fit 


vs. Log(L ) 


coefficient 


objects 


a chance correlation 


slope 


Log(of) 


0.12 


107 


1. 


— 


Log(R OM ,) 


0.90 


60 


lO" 20 


0.3 


Log(M mv ) 


0.96 


46 


lO" 22 


1.3 


Log(nioooAt/) 


0.72 


56 


10" 8 


0.3 


Log(<n>) 


0.25 


39 


0.4 


— 


Log(nio^r) 


-0.26 


18 


0.5 


— 


Log(#) 


0.88 


60 


lO" 19 


-6xl0- 3 



Table 6. Linear correlation coefficients and probability for a chance correlation between the pair of parameters plotted on Fig. [T] and 

m 
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Fig. A.15. Observed brightness profiles at 450 yum (upper left panel), 850 jim (upper right panel) and the SED (lower centred panel) 
toward Serpens FIRS 1. The solid lines report the best-fit model (Table |4]). The dashed lines represent the beam pattern of the 
telescope at 450 jim and 850 jim. The vertical dotted lines show the ou ter radius of the best-fit model. The lower limits from Plateau 
de Bure interferometric fluxes at 1.3 mm and 3 mm (JFuente et al.ll2007l) are also reported on the SED as lower limit. 
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Fig. A. 16. Serpens FIRS 1 dust temperature (left panel) and H2 density (right panel) profiles from the best-fit model. The dashed line 
represents the difference between the dust temperature and the gas temperature derived using X(H20)/„ = lxlO -5 and X(H20) 0Ut = 
lxlO -7 (scale on right axis). 
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